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ABSTRACT 

The Degree Angular Scale Interferometer (DASI) has measured the power spectrum of the Cosmic Microwave 
Background anisotropy over the range of spherical harmonic multipoles 100 < I < 900. We compare this data, 
in combination with the COBE-DMR results, to a seven dimensional grid of adiabatic CDM models. Adopting 
the priors h > 0.45 and 0.0 < r^. < 0.4, we find that the total density of the Universe fl,o, = 1.04 ±0.06, and the 
spectral index of the initial scalar fluctuations = LOlilio g^, in accordance with the predictions of inflationary 
theory. In addition we find that the physical density of baryons flt/h^ = 0.022;!;Qgy3, and the physical density of 
cold dark matter flcdmh^ = 0.14 ± 0.04. This value of flhh^ is consistent with that derived from measurements of 
the primeval deuterium abundance combined with big bang nucleosynthesis theory. Using the result of the HST 
Key Project, h = 0.72 ± 0.08, we find that Q,ot = 1 .00 ± 0.04, the matter density Q,„ = 0.40 ±0.15, and the vacuum 
energy density J^a = 0.60 ±0.15. (All 68% confidence limits.) 



Subject headings: CMB, anisotropy, cosmology 

1. INTRODUCTION 

The angular power spectrum of the cosmic microwave back- 
ground (CMB) has much to teach us about the nature of the Uni- 
verse in which we live (Hu, Sugiyama, & Silk 1997). Measure- 
ments are improving rapidly (de Bernardis et al. 2000; Hanany 
et al. 2000; Padin et al. 2001), and for a wide variety of the- 
oretical scenarios the predicted spectra can be calculated accu- 
rately (Zaldarriaga & Seljak 2000). Comparison of the data and 
models allows quantitative constraints to be placed on the pa- 
rameters of the Universe in which we find ourselves, and is the 
subject of this paper 

The Degree Angular Scale Interferometer (DASI), along with 
its sister instrument the CBI (Padin et al. 2001) and the 
VSA (Jones 1996), is one of several new compact interferom- 
eters specifically designed for observations of the CMB. This 
paper is the third in a set of three. Paper I (Leitch et al. 2001) 
gives a detailed description of the instrument, observations, and 
data calibration. Paper II (Halverson et al. 2001) focuses on the 
extraction of the angular power spectrum from the calibrated 
interferometric data, and provides band-power estimates of the 
angular power spectrum of the CMB. In this paper we com- 
bine the low-Z measurements made by the COBE-DMR instru- 
ment (Bennett et al. 1996) with the new measurements from 
DASI to constrain the parameters of cosmological models. 

The layout of this paper is as follows. The considerations 
which drive our selection of the model and parameter space to 
probe are detailed in §2. In §3 we review the method used to 
compare band power data to theoretical spectra. The results of 
this comparison are described in §4, and in §5 we draw some 
conclusions. 



2. MODELS, parameters AND MODEL GRID 
CONSIDERATIONS 



Following the discovery of the CMB, and the realization that 
the Universe went through a hot plasma epoch, it was proposed 
that adiabatic density perturbations in that plasma would lead 
to acoustic oscillations (Peebles & Yu 1970), and a series of 
harmonic peaks in the angular power spectrum (Doroshkevich, 
Zeldovich, & Sunyaev 1978). It was assumed that the initial 
fluctuations were scale invariant only because this is the sim- 
plest possible case. It was not until later that (New) Inflation 
was proposed (Guth & Pi 1982; Bardeen, Steinhardt, & Turner 
1983; Hawking 1982; Starobinsky 1982) — an elegant cosmo- 
genic mechanism which naturally produces such conditions. 
The simplest versions of this theory also make the firm pre- 
diction that the Universe is exactly flat, i.e., has zero net spatial 
curvature. 

Although Inflation sets the stage at the beginning of the 
plasma epoch it has nothing to say about the contents of the 
Universe. Over the last several decades a wealth of evidence has 
accumulated for the existence of some form of gravitating mat- 
ter which does not interact with ordinary baryonic material; the 
so-called cold dark matter (CDM). Conflicting theoretical ex- 
pectations and experimental measurements led to the proposal 
that a third component is present — an intrinsic vacuum energy. 
This three component model is the scenario we have chosen to 
consider 

The density required to produce zero net spatial curvature is 
referred to as the critical density pc = 3Hq/(^ttG), where Hq 
is the Hubble constant {Hq = 100/j km s"' Mpc"'). It is con- 
venient to measure the present day density of a component of 
the Universe in units of pc, denoting this fi,: the density of 
baryonic matter ilt,, the density of CDM flcdm, and the equiva- 
lent density in vacuum energy JIa. Thus the density of matter 
is given by = il/, + il.cdm and the total density is given by 
f2,o, = rim + r^A. Since fi,- = pi/pc note that ilih^ is a physical 
density independent of the value of the Hubble constant. 
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To generate theoretical CMB anisotropy power spectra we 
have used version 3.2 of the freely available CMBFAST pro- 
gram (Zaldarriaga & Seljak 2000). This is the most widely used 
code of its type, and versions 3 and greater can deal with open, 
flat and closed universes. CMBFAST calculates how the initial 
power spectrum of density perturbations is modulated through 
the acoustic oscillations during the plasma phase, by the ef- 
fects of recombination, and by reionization as the CMB photons 
stream through the Universe to the present. The program sets 
up transfer functions, taking as input il/,, ^laim and Qa, as well 
as Ho, the optical depth due to reionization (tc), and some other 
parameters. It can then translate initial perturbation spectra into 
the mean angular power spectra of the CMB anisotropy which 
would be observed in such a universe today. Inflation predicts 
that the initial spectrum is a simple power law with slope close 
to, but not exactly, unity. 

In any given comparison of data to a multidimensional model 
we may have external information about the values of some or 
all of the parameters. This may come from theoretical preju- 
dice, or from other experimental results. It may also be that the 
data set in hand is unable to simultaneously constrain aU of the 
possible model parameters to the precision which we desire. In 
such cases we can choose to invoke our external knowledge, 
and fold additional information about the preferred parameter 
values into the likelihood distribution. Often this occurs be- 
cause a parameter which could potentially be free is fixed at 
a specific value (an implicit prior). Or we may multiply the 
likelihood distribution by some function which expresses the 
values of the parameters which we prefer (an exphcit prior). 
The choice of measure, e.g., whether a variable is taken to be 
linear or logarithmic, is also an implicit prior — although such 
distinctions become less important as a variable becomes in- 
creasingly well constrained. 

There is no a priori reason to suppose that the marginal like- 
lihoods of the cosmological parameters are Gaussian. Thus the 
curvature of the likelihood surface at the peak does not fully 
characterize the distribution. To set accurate constraints it is 
therefore necessary to explore the complete likelihood space by 
testing the data against a large grid of models. If the introduc- 
tion of priors is to be avoided the grid must be expanded until 
one can be confident that it encompasses essentially aU of the 
total likelihood. 

Unfortunately it turns out that even within the paradigm of 
adiabatic CDM models the anisotropy power spectrum of the 
CMB does not fully constrain the parameters of the Universe. 
For example, it is well known that il„, and J^a are highly degen- 
erate. It is always necessary to invoke external information in 
any constraint-setting analysis. The clear articulation of these 
priors, both implicit and explicit, is critically important, as has 
previously been noted (Jaffe et al. 2001). 

We have chosen to consider a seven dimensional model 
space. The parameters which we include are the physical densi- 
ties of baryonic matter (flbh^ = uit) and CDM (Slcdmh^ = i^cdm), 
as well as the spectral slope of the initial scalar fluctuations 
(«j), and the overall normalization of the power spectrum as 
measured by the amplitude at the tenth multipole (Cio = 1(1 + 
l)Cio/27r). Noting that the degeneracy in the {i}m,^A) plane 
is along a fine of constant total density we opt to rotate the 
basis vectors by 45° and grid over the sum and difference of 
these parameters: Q^ + Qa = and Q^-ilA = ^a- This 
reduces the size of the model grid required to box in the region 
of significant likelihood. The seventh parameter is the opti- 



cal depth due to reionization (tc). Note that for each point in 
(Q A , fifot , i^b , ^cdm) space there is an implied value of the Hub- 
ble constant (h = a/ 2(wi + cucdm) / i^ioi + so we can calcu- 
late the (Uh,ncdm,^A,h) values for input to CMBFAST. 

We assume, as is the theoretical prejudice, that the contribu- 
tion of tensor mode perturbations is very small as compared to 
scalar, and ignore their effect (Lyth 1997). Tensor modes pri- 
marily contribute power at low I numbers, so if a large fraction 
of the power seen by DMR were caused by this effect the scalar 
spectrum would need to be strongly tilted up to provide the 
observed power at smaller angular scales. However we know 
that Hs cannot be » 1 as this would conflict with results from 
large scale structure studies. Our constraints should, however, 
be taken with an understanding of our assumption regarding 
tensor modes. 

In principle some of the dark matter could be in the form of 
relativistic neutrinos (hot as opposed to cold dark matter). How- 
ever the change that this would make to the CMB power spec- 
trum is negligible compared to the uncertainties of the DASI 
data (Dodelson, Gates, & Stebbins 1996). We therefore assume 
that all the dark matter is cold, and set Q,,, = 0, although it should 
then be understood that the flcdmh^ value we find may, in prin- 
ciple, contain some hot dark matter. 

Papers fitting tiie BOOMERANG-98 and MAXIMA- 1 band- 
power data (Balbi et al. 2000; Lange et al. 2001; Jaffe et al. 
2001) considered seven dimensional grids rather similar to our 
own. Other studies have examined model grids with as many 
as 11 dimensions (Tegmark, Zaldarriaga, & Hamilton 2001), 
including an explicit density in neutrinos and tensor mode per- 
turbations, generally finding both effects to be small. 

Taking the philosophy that simplicity is a virtue, we have 
not made use of /-space or fc-space splitting to accelerate the 
calculation of the model grid (Tegmark et al. 2001). In ad- 
dition we have generated a simple regular grid, rather than at- 
tempting to concentrate the coverage in the maximum likeli- 
hood region. Finally we have not treated the normalization pa- 
rameter Cio as continuous, preferring to exphcitly grid over this 
parameter as well. This is computationally somewhat slower, 
but makes the marginalization simpler, and involves no assump- 
tion about the form of the variation of versus this parameter. 
Using the notation lower-edge:step-value:upper-edge (number- 
of-values) our grid is as follows: = -1-0 : 0.2 : 3.4 (23), 
ntot = 0.7 : 0.05 : 1.3 (13), Qbh^ = 0.0100 : 0.0025 : 0.0400 
(13), n^dmh^ = 0.00 : 0.05 : 0.5 (11), t, = 0.0 : 0.1 : 0.4 (5), 

= 0.75 : 0.05 : 1.25 (11), and Cio = 300 : 50 : 1300 (21). Ex- 
cluding the small, physically unreasonable corner of parame- 
ter space where flm < 0, we make 205,205 runs of CMBFAST 
generating 205,205 x 11 = 2,257,255 theoretical spectra and 
calculating 2,257,255 x 21 = 47,402,355 values of against 
the data. 

For theoretical and phenomenological discussions of how the 
various peak amplitudes and spacings of the power spectrum 
are related to the model parameters see Hu et al. (1997) and 
Hu et al. (2001). In this paper we choose to compare data and 
models without explicitly considering such connections. 

3. COMPARISON OF DATA AND MODELS 

Consider a set of observed band-powers Vj in units of /xK^, 
together with their covariance matrix Pjj. If the overall frac- 
tional calibration uncertainty of the experiment is s we can add 
this to the covariance matrix as follows: 

Nij = Pij + s^ViVj. (1) 
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For the purposes of the present analysis we assume s = 0.08 
which includes both temperature scale and beam uncertainties 
(the fractional error on the band-powers is 8% corresponding to 
4% in temperature units — see Paper II). 

Now consider a model power spectrum C/. The expectation 
value of the data given the model is obtained through the "band- 
power" window function /I (Knox 1999), 



Bond et al. 2000), concatenating the and x, vectors and 
forming a block diagonal covariance matrix. Note that while 
the effect of the transformation described above is modest for 
the DASI points, it is very important for those from DMR (due 
to the large sample variance at the lower /'s). 



(2) 



The band-power window functions W^/l are calculated from 
the band-power Fisher matrix, F, and the Fisher matrix of 
the bands, fo„ sub-divided into individual multipole moments, 

?=Er').E^«" (3) 

(adapted from Knox 1999, for Fisher matrices with significant 
off-diagonal elements). The sum of each row of the array / / 
is unity, so Equation (2) simply represents a set of weighted 
means. Note that any experiment with less than full sky cov- 
erage will always have non top-hat band-power window func- 
tions. In practice, we calculate Equation (3) by subdividing 
each band into four sub-bands, and interpolate the results. The 
functions for the DASI band-powers are plotted in Figure 1, and 
are available at our website'. In practice the effect of using the 
correct window function, versus simply choosing the C/ at the 
nominal band center, is extremely modest. 

The uncertainties of the "D, are non-Gaussian so it would not 
be correct to calculate at this point. However it is possible 
to make a transformation such that the uncertainties become 
Gaussian to a very good approximation (Bond, Jaffe, & Knox 
2000). An additional set of quantities need to be calculated 
from the data which represent the component of the total uncer- 
tainty which is due to instrument noise. We can then transform 
each of the variables as follows, 

Pf = ln(P,+x,) (4) 

T^ = ln(%+xd (5) 

Mfj=N;/iVi+Xi)iVj+xj), (6) 

and calculate as usual, 

X^ = (Vf-%^)Mfj(V^-Tf). (7) 

The inverse covariance matrix elements M? will be approx- 
imately independent of the bandpowers D,. This is true even 
with the added calibration uncertainty term in Equation (1), un- 
der the assumption that the band-power uncertainty is sample 
variance dominated, i.e., Vi/xi ^ 1 (as is the case with iilmost 
aU the DASI band-powers), or that the fractional caUbration un- 
certainty is small compared to the total uncertainty in the band- 
power, s-^ <^ Nii/Vf, which is the case for DMR. Use of this 
transformation is very important as it allows us to use x^> and 
therefore not only to find the best fitting model, but to determine 
an absolute goodness of fit. 

The ability of smaller angular scale (/ > 100) CMB data to 
set constraints on model parameters is much improved when 
the large angular scale (I < 25) information from the DMR in- 
strument is included. We use the DASI bandpowers described 
and tabulated in Paper 11 together with the 24 DMR band- 
powers provided in the RADPACK distribution (Knox 2000; 

' http : //astro . uchicago . edu/dasi 

^ In fact since there is not sufficient freedom in the model that an arbitrary set 
than 7. For example assuming 4 shifts us to the 56% point of the cdf. 



4. RESULTS 

Figure 2 shows the DASI band-powers, together with the 
DMR data condensed to a single point for display. The x^ 
of the best fit model which falls on our grid is 29.5 for the 
9 DASI plus 24 DMR band-powers. Assuming a full 7 de- 
grees of freedom are lost to the fit, this is at the 71% point 
of the cumulative distribution function (cdf)^. The param- 
eters of this model are (Q„,i}\,il,hh^,il,cdmh^,Tc,ns,CiQ) = 
(0.725,0.325,0.0200,0.15,0.0,0.95,800), equivalent to 
{VLb,Q.cdmMK,Tc,ns,h) = (0.09,0.64,0.33,0,0.95,0.48). How- 
ever, no particular importance should be ascribed to these — 
the concordance model (Ostriker & Steinhardt 1995; Krauss & 
Turner 1995) which is shown has a x^ of 30.8 (76%) and is 
also rather a good fit. The message of Figure 2 is simply that 
there are models within the grid which fit acceptably well, and 
that we are therefore justified in proceeding to marginalized 
parameter constraints. We convert x^ to likelihood, C = 1'^. 

The extreme degeneracy of CMB data in the (fi^, Qa) plane 
has already been mentioned. This inability to choose between 
models with the same ^tot is in fact weakly broken at very low-/ 
numbers by the Sachs-Wolfe effect (Efstathiou & Bond 1999). 
The likelihood contours diverge from the ^tot = 1 line as 
becomes ^ 1, and the allowed region broadens. Consequently 
the marginal likelihood curve of Cl,ot acquires a high side tail as 
models with progressively greater Ha are included. These high 
JIa models have very low values of h, and are known to be in- 
valid from a wide range of non-CMB data. This being the case 
it is clearly not sensible to allow them to influence our results. 

We are therefore prompted to introduce additional external 
information. We could simply restrict fi„, and n\ to some 
"reasonable" range; for example requiring J^a > and il„, < 
1. Instead we choose to introduce a prior on h since this 
strongly breaks the degeneracy and is a quantity which is be- 
lieved to have been measured to 10% precision (Freedman et al. 
2000). We use two h priors; a weak h > 0.45, and a strong 
/i = 0.72 ± 0.08 assuming a Gaussian distribution. For the weak 
prior adding the additional restriction h < 0.90 has almost no 
effect as the excluded models already have very small likeli- 
hoods. To apply the prior we simply calculate the h value at 
each grid point, assign it the relevant likelihood, and multiply 
the two grids together. 

Figure 3 shows how the marginal likelihood distributions of 
the model parameters change as we move from the implicit 
prior of (JIa < 3.4,rifof < 1.3), to weak, and then strong prior 
cases. Note that most of the curves fall to a small fraction of 
their peak value before the edge of the grid is reached. For all 
parameters where this is not the case one must introduce ex- 
ternal priors such that it becomes so, and/or acknowledge the 
implicit top-hat prior which the range of that grid parameter 
represents — only then can the constraint on any of the param- 
eters be accepted. AU such priors must then be quoted with 
the constraints. In fact we see that i^lbh^ ,Hcdmh^ ,ns,CiQ) are 

of 7 band-powers could be fit exactly the effective loss of degrees of freedom is less 
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Fig. 2. — The DASI first-season angular power spectrum in nine bands (closed circles). The DMR information is shown compressed to the single lowest / 
point. The solid (red) line is the best fitting model which falls on our grid, while the dashed (green) shows the concordance model (Qt^^cdnn^A^TcnsJi) = 
(0.05,0.35,0.60,0, 1.00,0.65). The error bars plotted here are strictly for illustrative purposes only. The calculation is made using the full covariance matrix, and 
after the transformation described in §3. Thus "Chi-by-eye" can be misleading. 
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Fig. 3. — Marginal likelihood distributions for each dimension of the model grid. The dotted (blue) and solid (black) lines show the distribution before and after 
the introduction of a weak prior on the Hubble parameter (h > 0.45). The dashed (red) lines shows the effect of the stronger prior h = 0.72 ± 0.08. In the Q^/i^ panel 
the BBN constraint is shown as a (green) shaded region. All curves are normalized to a peak height of unity, and are spline interpolations of the actual model grid 
values shown by the points. 



almost completely unaffected by the choice of prior on h — 
this indicates that correlations between these parameters and 
(r^A , ^lor) are modest, and is a fortunate result. 

We turn now to r^. which, as is evident in Figure 3, is a 
very poorly constrained parameter The effect of reionization 
is to suppress power at small angular scales, and tilt the spec- 
trum down. However this can be compensated by adjusting 
tis upwards making these two parameters largely degenerate. 
Worse still, unlike the h prior discussed above, we have very 
weak experimental guidance as to the value of — we know 
only that reionization occurred at z > 6, roughly equivalent to 
Tc > 0.03. From theoretical ideas regarding early structure for- 
mation, and energy emission, it seems essentially impossible 
that Tc > 0.4 (see Haiman & Knox 1999, for a recent review of 
our knowledge regarding reionization). The theoretical preju- 
dice for « 1 is strong, but since this is a fundamental param- 
eter of the cosmology which we are trying to measure, we are 
very reluctant to place a prior on it. 

We have opted to accept the top-hat prior implicit in our 
model grid, i.e., that the likelihood of Tc falling between 0.0 
and 0.4 is uniform. In Table 1 we list the spline interpolated 
median, Icr and 2a points of the integral constraint curves. The 
modal (maximum likelihood) value is also given. All of the 
constraints quoted in this paper are 68% confidence intervals 
about the median. Although one can argue that the mode is 
perhaps more natural, in practice it makes little difference. 

Referring to Figure 3 we see that the Ha constraint is almost 
wholly driven by the prior on h; for this reason we have not 
included it in table 1. However, if one is prepared to accept 
the strong prior h = 0.72 ± 0.08, then our data indicates that 
n„„ = 1 .00 ± 0.04, n„, = 0.40 ± 0. 15 and = 0.60 ±0.15. 

Figure 4 illustrates the effect of setting = 0.0 (which the 
data weakly prefers). As already mentioned the primary degen- 



eracy is with n, which shifts to = 0.97;|;Q{j4. 

The selection of the particular set of nine bandpowers which 
we have used in this analysis is quite arbitrary. We have tested 
increasing the number of bins — as expected the variance and 
covariance of the points increases to compensate, and the con- 
straint curves remain the same. Shifting the bins in I also 
leaves the results unchanged — for instance the alternate bin- 
ning shown in Paper II leads to results which are indistinguish- 
able from those presented here. 

5. CONCLUSION 

We have compared the DASIh-DMR data to adiabatic CDM 
models with initial power law perturbation spectra. The best fit- 
ting model has an acceptable value, indicating that for data 
of the present quality models within this class are a plausible 
representation of the underlying physics. Adopting the con- 
servative priors h > 0.45 and 0.0 < < 0.4, we find n,ot = 
1 .04±0.06 and n, = 1 .Olilio ge. consistent with the predictions of 
Inflation. Moving to more aggressive priors on h and Tc tight- 
ens the constraints on il,,,, and respectively but they remain 
consistent with the theory. 

We find that ^bh^ = 0.022!j;{jg4 and n^dmh^ = 0.14 ±0.04, 
adding to the already very strong evidence for non-baryonic 
dark matter These constraints are only weakly affected by the 
choice of h and priors. Setting a strong h prior breaks the 
(n„,,il,A) degeneracy such that we constrain fi,,, = 0.40 ±0.15 
and = 0.60 ± 0.15 — consistent with other recent results. 

The current best value for fthh^ derived from the well de- 
veloped theory of big bang nucleosynthesis (BBN), combined 
with measurements of the primeval deuterium abundance, is 
nbh^ = 0.020 ±0.002 (Buries, NoUett, & Turner 2001, 95% 
confidence). The of the difference between this and our 
own value is at the 42% point of the cdf (assuming Gaus- 
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Table 1 

Parameter constraints from DASI+DMR data 



Parameter 


2.5% 


16% 


50% 


84% 


97.5% 


mode 




0.927 


0.986 


1.044 


1.103 


1.150 


1.047 




0.0156 


0.0187 


0.0220 


0.0255 


0.0292 


0.0220 




0.075 


0.100 


0.137 


0.175 


0.225 


0.135 




0.901 


0.949 


1.010 


1.092 


1.166 


0.993 




558 


642 


741 


852 


973 


728 



These constraints are from a seven dimensional grid, assuming the weak prior h > 0.45, and 0.0 < < 0.4. The indicated points on 
the cumulative distribution function are given, as well as the maximum likelihood value. 




Fig. 4. — Marginal likelihood distributions when varying the prior on Tc. All curves assume the weak prior h > 0.45. The solid (black) are the same as in Figure 3 
and assume 0.0 < Tc < 0.4, while the dashed (red) set Tc = 0.0. 



7 



sian errors on both); the values are hence consistent. Previous 
CMB analyses have seen httle power in the second peak region, 
and have determined Q^h^ values higher than, and inconsistent 
with, BBN at the w 3a level (Lange et al. 2001; Jaffe et al. 
2001). 

We would like to thank the staff of Argonne National Labo- 
ratory for allowing us to use the "Chiba City" computer cluster 
to generate the model grid used in this paper. Lloyd Knox, 



Wayne Hu and Mike Turner are thanked for useful conversa- 
tions. Most importantly we thank U. Seljak & M. Zaldarriaga 
for making CMBFAST publicly available. We would also like 
to thank the CBI team at Caltech for assistance with the design 
and implementation of the DASI instrument. This research is 
supported by the National Science Foundation under a cooper- 
ative agreement (NSF OPP 89-20223) with CARA, a National 
Science Foundation Science and Technology Center. 



REFERENCES 



Balbi, A., Ade, R, Bock, J., et al., 2000, ApJ, 545, LI, astro-ph/0005124 
Bardeen, J. M., Steinhardt, R J., & Turner, M. S. 1983, Rhys. Rev. D, 28, 679 
Bennett, C. L., Banday, A. J., Gorski, K. M., et al., 1996, ApJ, 464, LI 
Bond, J. R., Jaffe, A. H., & Knox, L. 2000, ApJ, 533, 19, astro-ph/9808264 
Buries, S., Nollett, K. M., & Turner, M. S. 2001, Phys. Rev. D, 63, 063512, 
astro-ph/0008495 

de Bernardis, R, Ade, R A. R., Bock, J. J., et al., 2000, Nature, 404, 955 
Dodelson, S., Gates, £., & Stebbins, A. 1996, ApJ, 467, 10. astro-ph/9509147 
Doroshkevich, A. G., Zeldovich, I. B., & Sunyaev, R. A. 1978, AZh, 55, 913 
Efstathiou, G. & Bond, J. R. 1999, MNRAS, 304, 75 
Freedman, W. L., Madore, B. R, Gibson, B. K., et al., 2000, ApJ, 553, 47, 

astro-ph/00 12376 
Guth, A. H. & Pi, S. . 1982, Physical Review Letters, 49, 1110 
Haiman, Z. & Knox, L. 1999, in Microwave Foregrounds, ed. A. de Oliveira- 

Costa & M. Tegmark, Vol. 181 (San Francisco: ASP), 227, astro-ph/99023 1 1 
Halverson, N. W., Leitch, E. M., Pryke, C., et al., 2001, ApJ, 568, 38, 

astro-ph/0104489 

Hanany, S., Ade, R, Balbi, A., et al., 2000, ApJ, 545, L5, astro-ph/0005123 
Hawkmg, S. W. 1982, Phys. Lett. B, 115, 295 
Hu, W., Fukugita, M., Zaldarriaga, M., et al., 2001, ApJ, 549, 669, 
astro-ph/0006436 

Hu, W., Sugiyama, N., & Silk, J. 1997, Nature, 386, 37, astro-ph/9604166 



Jaffe, A. H., Ade, R A. R., Balbi, A., et al., 2001, Phys. Rev. Lett., 86, 3475, 
astro-ph/0007333 

Jones, M. E. 1996, in Moriond Astrophysics Meetings, Vol. XVI, Microwave 
Background Anistropies, ed. B. G. J. V. FR. Bouchet, R. Gispert 
(Gif-sur-Yvette: Editions Frontieres), 161, ISBN: 3863322087 
Knox, L. 1999, Phys. Rev. D, 60, 103516, astro-ph/9902046 
Knox, L. 2000, http://bubba.ucdavis.edu/ knox/radpack.html 
Krauss, L. & Turner, M. S. 1995, Gen. Rel. Grav 27, 1137, astro-ph/9504003 
Lange, A. E., Ade, R A. R., Bock, J. J., et al., 2001, Phys. Rev D, 63, 042001, 
astro-ph/0005004 

Leitch, E. M., Pryke, C, Halverson, N. W., et al., 2001, ApJ, 568, 28, 

astro-ph/0104488 
Lyth, D. H. 1997, Phys. Rev Lett., 78, 1861 

Ostriker, J. R & Steinhardt, R J. 1995, Nature, 377, 600+, astro-ph/9505066 
Padin, S., Cartwright, J. K., Mason, B. S., et al., 2001, ApJ, 549, LI, 

astro-ph/00I2211 
Peebles, R J. E. & Yu, J. T. 1970, ApJ, 162, 815 
Starobinsky, A. A. 1982, Phys. Lett., B117, 175 

Tegmark, M., Zaldarriaga, M., & Hamilton, A. 2001, Phys. Rev. D, 63, 

043007, astro-ph/0008167 
Zaldarriaga, M. & Seljak, U. 2000, ApJS, 129, 431, astro-ph/9911219 



